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Abstract 

In the past two decades high-harmonic generation (HHG) has become a key process 
in ultra-fast science due to the extremely short time-structure of the underlying 
electron dynamics being imprinted in the emitted harmonic light bursts. After dis- 
cussing the fundamental physical picture of HHG including continuum-continuum 
transitions, we describe the experimental progress rendering HHG to the unique 
source of attosecond pulses. The development of bright photon sources with zep- 
tosecond pulse duration and keV photon energy is underway. In this article we 
describe several approaches pointed toward this aim and beyond. As the main bar- 
riers for multi-keV HHG, phase-matching and relativistic drift are discussed. Routes 
to overcome these problems are pointed out as well as schemes to control the HHG 
process via alterations of the driving fields. Finally, we report on how the investi- 
gation of fundamental physical processes benefits from the continuous development 
of HHG sources. 
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1 Introduction 



The photon can be seen in many regards as an engine of technological progress 



for the 21st century (National Research Council and Committee on Optical 



Science and Engineering, 1998) — a success story that began with the advent 



of lasers about 50 years ago (Maiman, 1960). Technology has continuously 



developed and lasers have become an important tool in industry and in almost 
every physics laboratory, independent of its field of research. Besides other 
parameters, the achievable laser intensities are a significant means to depict 



the progress in laser technology (Mourou et al. 2006). The first lasers in the 



1960s exhibited intensities well below 10 10 W/cm 2 . An enormous interest of the 
scientific community in this tool and vivid research led to rapid progress of the 
highest laser intensities that can be accomplished. Employing chirped-pulse 
amplification (Strickland and Mourou, 1985[ ) renders it nowadays possible to 
create laser pulses with peak intensities of 10 22 W/cm 2 (Yanovsky et al. , 2008 ). 
Intensities of more than 10 24 W/cm 2 are envisaged for the upcoming Extreme 
Light Infrastructure (ELI) within the next decade, with probable impact on 
many areas in physics such as laser plasma acceleration, attoscience, warm 
dense matter, laboratory astrophysics, laser fusion ignition, nuclear and high- 
energy physics with lasers. Besides the high peak intensities, pulse durations 



down to only a single cycle have also been achieved ( 


Nisoli et al. 


1996; 


Nisoli 


et al. 


1997 


Schenkel et al. 


2003 


Wirth et al. 


2011), the carrier envelope 


phase can be stabilized ( 


Baltuska et al. 


2002 
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3altuska et al. 




2003 


), lasers 



have become more reliable and are available at a variety of wavelengths. 

With the invention of the laser, a bright coherent light source became available 
that had, among others, huge impact on atomic physics opening the field of the 
laser spectroscopy and led to a renaissance in the field of optics with the birth 



of non-linear optics (Franken et al. , 1961 Boyd, 2008). A substantial advance- 



ment in non-linear optics and atomic physics happened when the laser field 
strengths became comparable with the electric field experienced by a bound 
electron which was achieved in the 1980s. At that point, the non-perturbative 
regime of non-linear optics was entered. Fields with such a strength are able to 
significantly change the electronic dynamics of atoms and molecules. In par- 
ticular, large fractions of the electronic wave function can be transferred to the 
continuum within very short times. A number of new effects were discovered 



at that time. First we mention above-threshold ionization (ATI) (Agostini 
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Fig. 1. (Color inline) Schematic of the three-step model for the HHG process. The 
blue surface is the superposition of the atomic binding potential and the electrostatic 
potential of the laser field. The bound and continuum part of the wave function are 
rainbow-colored. A tiny part of the wave function (continuum part) is just freed 
and subsequently driven away. When the laser field reverses its sign, the continuum 
part is stopped and accelerated back towards the core. At recollision both parts of 
the wave function interfere and give rise to coherent photon emission denoted by 
the blue wiggled line. 



et al. , 1979) happening when electrons ionized by an infrared (IR) laser field 



were detected with a large kinetic energy well above the ionization threshold. 
Moreover, recollisions of those electrons with the core were observed, leading 
either to the knock-out of a second electron, the non-sequential double ion- 
ization (L'Huillier et al. 1983 Fittinghoff et al. 1992), or to recombination 



along with the emission of extreme-ultra- violet (XUV) light, called high-order 
harmonic generation (HHG) - the main process under consideration in this 



review. 



HHG was discovered in 1987 (McPherson et al. , 1987; Ferray et al. , 1988): 



a rare gas illuminated by a laser emitted photons with an energy of several 
odd multiples of the laser frequency. Remarkably, the emission spectrum ex- 
hibits a plateau-like structure extending far until a sudden cutoff rather than a 
continuous exponential decrease typical for perturbative non-linear processes. 
This rapidly led to intensive discussions about the physical origin of the phe- 
nomenon which eventually brought forth the three-step model (see Fig. [IJ 



discovered in 1993 (Corkum 1993; Schafer et al. 1993). i) the electronic wave 



function of an atom is partially freed by a strong laser field, ii) then the ionized 
fraction is subsequently driven away in the continuum by the laser field, hi) 
finally, the wave packet is accelerated back to the ionic remnant, interfering 
with the bound part of the wave function giving rise to a strong, coherent, 
high-frequency response that can lead to the emission of a HHG photon along 
with the recombination of the electron into the bound state. 



HHG is fascinating both from being a fundamental example of nonperturba- 
tive laser-atom interaction dynamics and from a technological point of view, 
allowing creation of light flashes with exceptional properties: the energy of an 
emitted HHG photon is the sum of the binding energy and the kinetic energy 
acquired in the laser field. The photon energy can be extremely large reaching 
several kilo-electronvolts (keV) ( |Seres et al. , 2006). In a world where techno 1- 
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ogy is currently miniaturizing, x-rays with their tiny wavelengths may become 
the light of the future and thus HHG potentially of particular importance. 
Another property of HHG is its coherence, which arises because the phase 



of harmonic emission is locked to the laser phase (Lewenstein et al. 1995), 
and under optimal conditions the harmonic yield scales quadratically with 
the number of contributing atoms. The last point indicates that the process 
cannot be viewed only from a single-atom perspective. It is crucial for the emis- 
sion that all atoms are phase-matched meaning that harmonics from different 
atoms have to emit with the same phase to allow for constructive interference 
between them. While the macroscopic aspect of HHG (phase-matching) is im- 
portant to get a sizable harmonic yield and to create XUV/x-ray sources, the 
single-atom aspect is not less essential. Thus, the harmonic spectrum contains 
structural signatures of the emitting atom or molecule which can be assessed 



by means of HHG spectroscopy ( 


Lein et al. 




2002 


Nalda et al. , 


2004 


Itatani 


et al. 


2004 


Lein 


2005 


Morishita et al. , 


2008 


)• 



HHG opened the door for attoscience (Krausz and Ivanov, 2009). In fact, the 



bandwidth of emitted harmonics is large enough to allow for generation of light 



2006 Goulielmakis et al. , 2008 Ko et al. , 2010). Attoseconds is the time scale 



bursts much shorter than the laser pulse duration. The harmonic pulses ( 


Nisoli 


and Sansone, 2009 Chang and Corkum 


2010 


Sansone et al. , 


201 1| Calegari 


et al. 


2011 ) can have durations down to the attosecond regime ( 


Sansone et al. 



of electron dynamics in atoms and ions. For a long time, no tools have been 
available that were fast enough to resolve such dynamics. For this reason, HHG 
has become a key process in atomic and molecular physics and a state-of-the- 
art coherent XUV light source being nowadays available in many research 
laboratories. With the emergence of extremely short pulses from HHG, an 
entirely new pathway in exploring and controlling the physics of atoms and 



ions ( 


Scrinzi et al. 


2006 


jKienberger et al. , 


2007 


Nisoli and Sansone 


2009 


Krausz and Ivanov 


2009 


Brif et al. , 


2010) was opened. 



The capability of the current HHG sources is limited to energies around a 
few hundred eV and to pulse durations of several tens of attoseconds. An 
advancement to shorter pulse durations and higher photon energies is highly 
desirable because it will allow not only to track the electron dynamics in atoms 
but also to time-resolve intra-nuclear dynamics. In this article we present 
selected approaches to reach these aims. For previous reviews on similar topics, 



we refer the reader to 


Protopapas et al. 


(199' 


D;S 


alieres et al. 


(1999); 


Brabec 


and Krausz (2000 


); 


Joachain et al. (200C 


)); 


Becker et al. 


(2002 


; Agostini and 


DiMauro (2004); 


Pfeifer et al. 


(2006c); 


Salamin et al. ( 


2006); 


Mourou et al. 


(2006 


); 


Scrinzi et al. 


(2006 


; Gaarde et al. 




2008 


; Win 


terfeldt et al. 


(2008); 


Nisoli and Sansone 


( 


2009) 


Krausz and Ivanov 


(2009) 


Popmintchev et al. 


(2010 


); 


Sansone et al. 


(2011 


>;< 


Dalegari et al. 


(201 


1). 



An alternative technique for HHG is the irradiation of solids with super-strong 



lasers (/ > 10 16 W/cm 2 ) leading to the formation of overdense plasmas and 



the light emission from its surface (] 


Jchters et al. 


1996 


Quere et al. 


2006 


Dromey et al. 


2009 


Nomura et al. 


2009 


). In our review we 


do not cover 


this approach and refer the reader to, e.g., 


Tsakiris et al. 


(2006 


); 


Thaury and 



Quere (2010). 



2 Fundamental concepts of HHG and attosecond pulses 



We begin with a brief presentation of the basic principles and the phenomenol- 
ogy of HHG along with an introduction to some particular theoretical con- 
cepts. This leads us to a discussion about the interference picture of HHG 
including continuum-continuum transitions. Further, we review milestone ex- 
periments that have led to the corresponding state-of-the-art experimental 
technology. 



2. 1 Lewenstein model and Phenomenology 



The first theoretical description of HHG spectra has been accomplished by Lewpn 



stein et al. (1994) on the basis of the strong-field approximation (SFA) ( |Keldysh 



1965 Faisal, 1973 Reiss, 1980). In this work, a quasi-classical trajectory-based 



picture of the HHG process has been established where the evolution of the 
wave function in the continuum, derived from quantum mechanics, can finally 
be described by an ensemble of classical trajectories. The theoretical model 
is developed as follows. The harmonic emission is evaluated via the electric 
dipole moment (in atomic units) [Jj 



d(*) = -<¥(f)|x|*(f)> 



(1) 



where the electron wave function \I/(t) is determined within the SFA adopting 



the single-active electron approximation (Kulander, 1988 Lewenstein et al 



1994 Paulus et al. 1994). The SFA is required for finding an analytic solu- 



tion of the initially bound wave function evolving in the strong laser field. 
The fundamental SFA assumptions are that the binding potential is dominant 
before ionization and the laser field after ionization. This way, the Fourier 
transformed dipole moment d as a function of the harmonic frequency wh was 
found to be 



In atomic units the electron mass is m e = 1, its charge is e = —1, Planck's 
constant is h~ = 1 and consequently the speed of light is found to be c = 137. 
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_ poo rt r 

d(u H ) = i dt dt' / d 3 p(0ok|p + A(*)/c) 

J — CO J -co J 

x (p + A(0 /c\xE{t) \fa) e -W(p,tt)-i 9 V-t)-wat) ; ( 2 ) 

where 

S(p,t,t') = l£ dr(p + A(r)/c) 2 (3) 

is the classical action of the laser field and I p the binding potential, |0o) 
the bound state and |p) the eigenstate of the momentum operator. A(t) and 
E(i) = —dtA(t)/c are the vector potential and electric field of the laser field, 
respectively. 

The integrand in Eq. ^ is highly oscillating because of the complex argument 
of the exponent and is numerically difficult to calculate exactly. Thus, one or 
several integrals can be carried out within the saddle-point approximation (see, 



e.g., Arfken and Weber (2005)) converting the quantum mechanical expression 
into a semi-classical expression. This means the integral is only evaluated 
around the stationary points of the phase S(p, t, t') — I p {t' — t) — u^t. These 
points are defined by the so-called saddle-point equations 



r dt"[ P +A(t")/c]=o 

[p + A(0/c] 2 /2 = -/ P 
[p + A(t)/c} 2 /2 + I p = io H 



(4) 

(5) 
(6) 



which correspond to energy conservation at ionization [Eq. ([5])] and recollision 
[Eq. (|6~j)] an d to the recollision condition of the classical electron [Eq. Q] (Lewen 



stein et al. , 1994). The stationary points or saddle points themselves represent 
the ionization and recollision time of the considered classical trajectory as well 
as the canonical momentum. It turns out that it is sufficient to sum over a 
small number of classically allowed trajectories for each energy to calculate the 



see, e.g., 

and for its relativistic extension Milosevic, D. B. et al. 



; Becker et al. 


(2002) 


(2000 


2002 


); 


Salamin 



et al. (2006); Klaiber et al. (2007a). Apart from the analytic approach above, 



the laser-atom-dynamics governed by the Schrodinger equation or relativistic 



wave equation can also be solved numerically (Muller, 1999 Bauer and Ko- 



val, 2006; Mocken and Keitel, 2008; Ruf et al. , 2009). Moreover, an approach 



based on the SFA but beyond the semi-classical picture has been developed 



recently (Plaja and Perez-Hernandez, 2007). 



The semi-classical picture has opened the perspective to intuitively under- 
stand a wide range of problems connected with HHG. Inspecting the classical 
trajectories, we briefly explain several principles of HHG. In Fig. [2] (a), we 
show different classical trajectories (colored lines) together with the laser field 
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(gray dashed line). Recollision is only possible for trajectories ionized in the 
quarter cycle after each maximum or minimum of the laser cycle. Trajectories 
starting before do no re-encounter the origin. All trajectories originating from 
the same half cycle recollide at different times and have different recollision 
energies. This also means that different harmonic wavelengths are emitted at 
different times and, thus, the emitted light has an intrinsic chirp, the so-called 



attochirp (Mairesse et al. 



2003). The recollision energies versus the recollision 



phase are shown in Fig. [2 (b). It can be seen that each energy is emitted twice 
per half cycle. These two branches of trajectories are termed long and short 
according to their excursion time in the continuum. As indicated in Fig. [2] (b), 
there is one trajectory with a maximum energy u c . This energy is called the 
cutoff energy and can be derived from the classical equations of motion yield- 



ing (Corkum, 1993) 



3.17 U n 



(7) 



where I v is the ionization energy and U p = t-% is the ponderomotive potential 
equal to the average quivering energy of an electron in a sinusoidal laser field 
with peak strength E and frequency u. 




Laser phase 



Recollision phase 



Fig. 2. (Color online) The gray dashed line shows the laser field in both parts of 
figure, (a) The solid colored lines are different classical trajectories. The trajectories 
start at different phases in the laser field. The recollision energy is encoded in their 
color. The energy increases from red to blue and the trajectory with the maximum 
energy (cutoff) is marked in black. In (b) we see the energies of the recolliding 
trajectories for different recollision times as solid line. 

A typical HHG spectrum is shown in Fig. [3] (a). The part of the spectrum 



below I p is called below-threshold harmonics (BTH). In the early years of 



laser physics (Franken et al. 1961) this branch of the spectrum has been 



modelled by taking the laser field perturbativly into account but it is under 



active discussion at the moment (Hu et al. , 2001; Yost et al. , 2009 



Power 



et al.[ |2010[ |Hostetter et al.[ |2010[ |Soifer et al.[ |2010} |Liu et al.[ |2011[ ). The 
long plateau in the spectrum ending at the cutoff energy uj c originates from a 
superposition of all classically allowed trajectories. The red line in (b) shows 
a state-of-the-art experimentally measured HHG spectrum ranging to keV 



photon energy, from Seres et al. (2005). Macroscopic effects are responsible 



for the deviation between (a) and (b). Moreover, the single harmonics cannot 
be resolved for such a large bandwidth. 



7 




Harmonic energy <o H Photon energy (keV) 



Fig. 3. (Color online) (a) Schematic of a typical HHG spectrum. The below-thresh- 
old harmonics (BTH) exhibit an exponential decrease in the first part followed by 
the plateau and a steep drop after the cutoff, (b) Experimentally measured HHG 
spectrum (red line, lowest line for high energies), transmittance function of the var- 
ious filters (green, highest line for high energies) and the theoretical single atom 



emission yield (gray, diagonal line). Figure reprinted with permission from Seres 



et al. (2005). Copyright 2005 by Nature (London). 



Unfortunately, HHG has a poor conversion efficiency. This is mainly because 
the free wave packet undergoes quantum spreading during the continuum mo- 
tion. When the wave packet is driven back, most parts of the wave function 
miss the core due to its large spatial extent after it has spread out. The im- 
pact of the spreading can be easily estimated. The dimensions of the bound 
wavefunction are typically on the order of la.u.. Th e transversal spread- 



ing velocity can be estimated by v x = v / 3^/(2/ p ) 1/4 (|Popov 



2004 



Ivanov 



et al. 1996). For a typical HHG experiment (E ~ 0.1 a. u., I p ~ la.u. and 



o;l = 0.057 a.u., where is the laser frequency), we find a wave packet di- 

30 a. u. at recollision. Thus, spreading in the two 



mension of x 



v± 



2tt 



perpendicular dimensions with respect to the laser polarization axis reduces 
the HHG yield by a factor of (30/1) 2 ~ 10 3 . 



Let us come back to the time structure of HHG. Around each maximum of the 
laser field, portions of the wave function tunnel out, recollide after a fraction 
of the laser period [see classical trajectories in Fig. [2] (a)] and lead to the 
formation of a light burst. In a periodic laser field the process repeats itself 
each half cycle of the driving field and the emitted field ]?h(wh) only differs in 
a phase-shift of n and the sign. Thereby it is assumed that two adjacent half 
cycles of the laser field are always identical apart from their sign and the atom 
possesses a central symmetry and thus the continuum dynamics are inverted 
after each half cycle leading to the difference in the sign of the dipole moment. 
A typical spectrometer being slow on the timescale of the laser period would 



thus measure a superposition of all these emitted pulses (Protopapas et al. 
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1996) 



L37T 



-Eh(^h) 



3 



•lVe 1 " 



'-JIT 



with significant contributions to the spectrum only at odd multiples of the laser 
frequency as in Fig. [3] even though the individual pulses exhibit a continuous 
spectrum. That is the physical reason for calling the effect high-harmonic 
generation. When violatating the former assumptions, even-order harmonics 



can also also be generated, e.g., in two-color laser fields (Kondo et al., 1996) 



via asymetric molecules (Gavrilenko and Oks, 2000; Kreibich et al. 2001) or 



driving pulses short enough that the periodicity is broken (Spielmann et al 



1997). Experimental progress meanwhile allows for the restriction of the HHG 
process to a half cycle only (see Sec. 2.4 for more details). In this case a 



continuous spectrum is measured. The duration of such a pulse (At) depends 
on two issues: first, the available spectral bandwidth (Ao>h) of the harmonics 
dictates the shortest possible pulse duration via the Fourier limit. It can be 
expressed by the time-bandwidth product 



Au YiAt = 3.6 as keV 



(9) 



under the assumption of a plateau-like spectrum. Second, the Fourier limit 
or bandwidth limit can only be reached when the pulse has no chirp. Thus, 
the intrinsic attochirp has to be compensated to compress the harmonic pulse 
toward its fundamental limit. 



2.2 Interference model of HHG 



One of the most fundamental properties of harmonic light flashes is their 
coherence, i.e., the phase of the harmonic light emitted from a single atom 
is locked to the laser phase and under suitable conditions even the phases of 
the harmonics from different atoms match. This particular property makes 
HHG a unique process with a variety of applications, for instance, rendering 
attosecond pulse generation possible. In order to understand why the light can 
be coherently emitted, the generation process of a single atom has to be seen 
from a quantum mechanical perspective: For a typical HHG scenario in the 
tunnel ionization or multiphoton regime (Protopapas et al. 1997), the initially 
bound wave function is partially ionized. The fraction of the wave function 
promoted to the continuum is then driven in the strong laser field and can 
eventually recollide. At that point, both, the ionized and the bound parts 
of the wave function interfere within the binding potential giving rise to a 
strong, coherent high-frequency dipole response that can lead to the emission 
of an HHG photon along with the recombination of the electron into the 



bound state (Pukhov et al. 2003 Itatani et al., 2004). We term the process 
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continuum-bound (CB) HHG because a continuum wave packet and a bound 
wave packet are involved. 



The interference picture of HHG can also be illustrated more quantitatively. In 



an electron-atom collision, the coherent part of the emitted radiation (Lappas 



et al. 1993 Burnett et al. 1992) can be calculated via the expectation value 



of the acceleration using the Ehrenfest theorem 

a(f) = -<*(f)|W|*(f)> (10) 

with the ionic potential V. Suppose the wave function would be a superposi- 
tion of a recolliding plane wave with momentum p and the bound state |</>o) 

■ p 2 

as in the usual HHG scenario: \^(t)) = ai|p)e -1 ^* + ci2|0o)e l/pt . The resulting 

p 2 

acceleration a(t) = — ala 2 {ei\ W\e2)e~ 1 ^ +Ip ' t + c.c. would both oscillate and 

2 

lead to emission at the expected frequency uu = \ + I p - In the case, the 
bound wave function has been depleted (a 2 = 0), the oscillations along with 



the coherent HHG emission are absent (Pukhov et al. , 2003). This situation 
occurs for laser fields strong enough to ionize the whole wave function before 
recollision. The intensity regime is called saturation or over-the-barrier ioniza- 
tion (OBI) regime and is reached when the barrier-suppression field strength 

j2 1 1 11 

-E'bsi = 4§ ( Augst et al.[ |1989[ ) is exceeded, where Z is the residual charge 



seen by the ionized electron. 

The former simple model based on eigenstates holds strictly true only for time- 
independent Hamiltonian but its picture remains true in atomic HHG with 
collision times typically short on the time scale of a laser period. Note that 
apart from coherent radiation considered in this review, a recolliding electronic 
wave packet can also cause emission of incoherent radiation due to spontaneous 
recombination. This spontaneous process also occurs when the electronic wave 
function has been fully depleted. The phase of the spontaneous recombination 
amplitude is random because of the arbitrary phase of the final state (ground 
state). However, in a macroscopic gas target this kind of emission yield scales 
only linearly with the number of atoms instead of the quadratic scaling typical 
for the coherent CB transition in the phase-matched case, in addition to it 
being non-directional as compared to the directed CB HHG light. 

Needless to say that further possibilities exist to create a time-dependent 



acceleration Eq. (10) leading to coherent photon emission. Especially relevant 
in HHG is the presence of two continuum wave-packet contributions within the 
range of the ionic potential. In this case, the wave function can be illustrated 

as 

|*i(*)> = ai|Pi)e-^ /a + a 2 \ P2 )e~^ 2t (11) 
leading to the acceleration a(t) = — a\a 2 {p\\ W\p 2 )e~ l( - P ^~ p ^^ 2t + c.c. that 
is beating and coherently emitting Bremsstrahlung at the difference energy 
ip 2 — Pi)/2 between both wave-packet fractions. This class of transitions is 
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termed continuum-continuum (CC) HHG and it has already been noticed in 



the early years of HHG (Lewenstein et al. , 1994). Although it was confirmed 



by several authors in the past (Becker et al. 1997 Kuchiev and Ostrovsky 



1999 Milosevic, D. B. 2001 Milosevic, D. B. and Ehlotzky 2003) that CC 



HHG has a magnitude of several orders of magnitude lower than CB HHG, 



interest in this transition has increased recently (Plaja and Perez- Hernandez 



2007 Perez-Hernandez and Plaja, 2007| Perez-Hernandez et al. , 2009| Kohler 



et al. , 2010) as ways have been found to separate it from the generally stronger 



CB HHG. 



On the bottom line of the former discussion, a unified picture for HHG includ- 
ing the CB and all CC transitions can be formulated: the interference of two 
parts of the same wave function with different energies in a potential causes 
coherent photon emission at their difference energy either due to recombina- 
tion radiation or due to Bremsstrahlung (Kohler et al. , 2010). For instance 



m 



the case of CB HHG, a recolliding electronic plane wave (energy €%) interferes 
with the bound wave packet (energy — I p ) leading to photon emission at the 
difference energy (ej + J p ). 



2.3 Continuum- continuum HHG 



As mentioned in the previous section, CC HHG has a several orders lower 
magnitude than CB HHG and, hence, its experimental discrimination is de- 
manding and has not been accomplished yet. The reason for this is the small- 
ness of the free-free transition cross-section (~ r%, with the electron classical 
radius r e ) with respect to the free-bound one (~ r|, with the Bohr radius r^). 
In the following we distinguish the three different CC transitions described in 
the literature and discuss new results which show the pathway to separate CC 
HHG from CB HHG. All CC transitions share that interference between two 
continuum wave packets within the range of the binding potential [see Fig. [2] 
is responsible for the emission of coherent Bremsstrahlung. The scenarios can 
be distinguished by the preceding dynamics of the involved wave packets. 



2.3.1 CC HHG with a single continuum wave packet 



The schematic in Fig. |4j (a) displays the recollision of a single wave packet 
with the ionic core. Since the ground state is not populated, no coherent (CB) 
emission takes place at the sum of the kinetic energy and the binding energy 
as argued before. However, due to the nonzero spectral width of the wave 
packet, the requirement of at least two simultaneously recolliding different 
momentum components is met [see Eq. (11)]. The recolliding electronic wave 
packet scatters at the ionic core which results in the transition between the 
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(a) 






Fig. 4. Schematic of the different CC transitions occurring in HHG. We concentrate 
on the coherent response. The thick black line is the potential, the thin gray lines 
are the different wave-packet portions. In (a) Bremsstrahlung is emitted when a 
single wave packet (CI) interacts with the core. In (b) the recolliding wave packet 
will interfere with the bound wave packet (B) and a just-emerging continuum wave 
packet (C2) leading to CB and CC harmonics, respectively. In (c), two continuum 
wave packets (C2 and CI) of different momentum recollide at the same time and 
interfere with each other emitting photons of exactly their kinetic-energy difference. 

two momentum states available in the electron wave packet along with photon 



emission forming coherent Bremsstrahlung (Protopapas et al. 1996 Watson 



et al. 1997 Milosevic, D. B. and Ehlotzky, 2003). The availability of ground- 
state population is irrelevant for this process and, therefore, it also occurs 
for laser intensities in the saturation regime where the whole electronic wave 
function is completely ionized. 



Very short radiation bursts down to the attosecond regime can even be emitted 



by this mechanism as calculations by Emelin et al. (2005, 2008) show. For 



this purpose, strong laser fields are required to ionize the electron within an 
extremely short time interval. In this case, the electron wave function is a 
single compact continuum wave packet after ionization. Later at recollision, the 
overlap between the wave function and the ion along with the emission process 
via the CC-transitions is confined to an attosecond time window under the 
precaution that the wave packet spread is within appropriate limits. According 
to the finite spectral width of the recolliding wave packet, the spectral range 
of the discussed CC transitions is limited and the photon energies are low 
compared to the kinetic energy at recollision and CB HHG. 



2.3.2 CC HHG with a Rydberg state 



Another kind of CC transition [shown in Fig. [4] (b)] happens when a rec- 
olliding wave packet interferes with a just-ionizing part of the bound wave 



function (Lewenstein et al. 1994 


Becker et al.[|1997| |Kuchiev and Ostrovsky 


1999 


Milosevic, D. B., 


2001| Mi 


osevic, D. B. and Ehlotzky, 2003| Plaja and 


Perez-Hernandez 2007 


Perez- Hernandez and Plaja[ 2007; Perez-Hernandez 


et al. 


2009) within the range of the ionic potential. As CB HHG, it only oc- 



curs in the saturation regime (with E < -E^si) where the atom is not depleted. 
This CC transition always accompanies CB HHG and is even indistinguish- 
able from the CB transition into an excited Rydberg state where the latter 
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has been populated from the initial ground state in the laser field. Only the 
energy of the CC transition is slightly smaller than that of the CB transition. 
Thus, the emission spectrum has similarities imprinted: the CC spectrum ex- 



hibits the same plateau-like structure as the CB spectrum, see Milosevic, D 



B. (2001). 



Due to the weakness of the CC spectrum with respect to CB HHG, it can 



be neglected in many cases (Kuchiev and Ostrovsky, 1999; Milosevic, D. B. 
2001). However, recent results (Plaja and Perez-Hernandez, 2007) indicate 



significant CC contributions in the multiphoton regime. In this work, the ratio 
between the CC and CB amplitude was estimated to be 1 — l/(47 2 ) with the 
Keldysh parameter 7 = JI P /(2U P ). 



2.3.3 CC HHG with two continuum wave packets 



In what follows, we consider two distinct continuum wave packets recolliding 
at the same time as sketched in Fig. H (c). In contrast to the transitions in 



Section |2.3.1| and |2.3.2[ both wave packets have been ionized in two different 
half cycles and have been evolved in the continuum for some time. Thus, both 
wave packets recollide with a kinetic energy that is nonzero and comparable 
with the ponderomotive potential of the laser field. The scenario can occur in 
multi-cycle laser pulses as sketched in Fig. [5] where classical electron trajecto- 
ries are shown (solid lines). Two portions of the wave function ionized at the 
starting times of the two trajectories will recollide simultaneously and result 
in CC emission at their kinetic energy difference. 
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Fig. 5. (Color online) Classical electron trajectories (solid lines) evolving in a mul- 
ti-cycle laser field (dashed line). 



In contrast to the CC transition in Section 2.3.2, this transition also occurs 
when the wave function is fully depleted at the recollision moment and, thus, 
can play a dominant role for HHG in the saturation regime. This regime has 



been studied by means of hydrogen atoms in Kohler et al. (2010) by a numer 



ical solution of the three-dimensional time-dependent Schrodinger equation 



using the Qprop of Bauer and Koval (2006). The laser pulse is shown in 
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Fig. [6] (a) and is chosen such that almost complete depletion of the ground 
state occurs on the leading edge of the pulse [see dashed line in Fig. M (a)]. To 




-10 



200 
time [a.u.] 



Fig. 6. (Color online) Time-frequency analysis of HHG showing the signature of 
CC wave-packet interference, a) The laser pulse used for the calculation (solid line, 
left axis) and ground-state population (dashed line, right axis), b) The windowed 
Fourier transform of the acceleration Eq. ( |10[ ) . The two dashed black lines represent 
the classically calculated kinetic energies of electrons returning to the ion and the 
solid red line represents their difference energy. Figure reprinted with permission 



from Kohler et al. (2010). Copyright 2010 by the American Physical Society. 



analyze the time-resolved frequency response of HHG, the windowed Fourier 
transform of the acceleration expectation value obtained from the TDSE cal- 
culation is calculated and displayed in Fig. [6)3. For comparison, the two dashed 
black lines in the figure display the classical recollision energies for trajectories 
starting from two different laser half cycles [first two peaks of the laser pulse] , 
respectively, which are in agreement with the traditional CB signal whereas 
the red line represents their difference. The CC transition is evidenced by the 
excellent agreement of the quantum-mechanical response with the red line. 
Interestingly, the CC component of the dipole response is the dominant con- 
tribution for several half cycles after t = 150 a.u.. This can be understood from 
the fact that depletion of the ground state occurs around that time. Then, co- 
herent HHG can only occur by the presence of the various parts of the wave 
function in the continuum. 



Moreover, the CC transition can be described within a strong-field approxi- 



mation model suitable for the OBI regime (Kohler et al.L 2010). The model is 



based on the evaluation of the acceleration a(t) = — {&(t)\ W\^(t)) ( |Gordon 
and Kartner , 2005b ) rather than the dipole moment to include the distortion 



of the recolliding waves by the Coulomb potential required for momentum con- 
servation. The saddle-point approximation is applied to the expression render- 
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Fig. 7. (Color online) Intensity dependence of ckcb and acc for the six quantum 
orbits at the harmonic frequency wh = 8 a.u. with respect to the peak pulse intensity 
which was chosen to be 10 16 W/cm 2 in the example of Fig. [6j The solid red lines 
represent the two CC contributions. The dashed black line depicts the values of ocb 
for the CB long and short contributions with a single return. The dotted black line 
denotes the CB contributions that arise from trajectories emitting at the second 
return. 



ing a computationally fast evaluation of the process possible. This analytical 
model also allows to extract the HHG dipole emission phase 



</>(!) = aj, (i = CB,CC) 



(12) 



being a function of the laser pulse peak intensity / (Lewenstein et al. , 1995) 



The results are shown in Fig. [7] for a photon energy of uh = 8 a.u. emitted 
between t = 160 a.u. and t = 210 a.u. where 6 different contributions exist [see 
Fig. ^jp] . The dashed and dotted lines represent «cb for the usual CB transi- 
tions with single and multiple return. The red solid line displays the intensity 
dependence of «cc f° r the CC transition. The result of Fig. [7] reveals a strik- 
ing difference for CB as compared to CC transitions: the sign of acc and c^cb 
differs for both types of transitions, consequently, the phase- matching condi- 
tions for macroscopic propagation will be different for these transitions. The 
different phase- matching behavior could allow for the discrimination of the CC 
harmonics from the CB harmonics after propagation through the medium. 



The intensity of CC HHG depends on the effective atomic or molecular (ionic) 
potential mediating the transition between the two continuum electron wave 
packets. This feature can be employed for qualitatively advancing tomographic 



molecular imaging (Itatani et al. , 2004): instead of probing the orbital shape 



of the active electron, the effective atomic or molecular potential could be 
assessed by investigation of CC spectra. 
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Fig. 8. (Color online) The situation of a coherent superposition of the emitted light 
(blue wiggled line) from different atoms (black circles) is shown. The HHG process 
of each atom in the Gaussian focus (red, gray area) is coherently triggered by the 
laser field which propagates to the right. 

2.4 Experimental advances 



Since the first discovery and experiments on high-order harmonic generation 



(HHG) in the late 1980s (McPherson et al. 1987 Ferray et al. 1988), a con- 



tinuous and still ongoing revolution in experimental technology has lead to 
progressively enhanced HHG and attosecond light sources as well as to the 
remarkable understanding and control over this nonlinear-optical light conver- 
sion process. The experimental setup is conceptually simple and involves only 
the focusing of a sufficiently high-energy ultrashort laser pulse (driver) into a 
gaseous medium. HHG takes places in each atom microscopically as discussed 
in the previous sections, and macroscopic phase matching of the individual 
emitters throughout the interaction region (which typically has lengths rang- 
ing from sub-mm to few cm and diameters in the few to hundreds of micron 
regime) ensures the generation of a bright coherent beam of harmonics after 
the conversion medium, typically copropagating with the generating funda- 
mental driver pulse. 

For phase matching, essentially, the vectorial sum of the individual n funda- 
mental photon wavevectors has to equal the wavevector of the nth-harmonic 
photon (i.e. photon momentum conservation) to ensure perfectly constructive 
interference of all microscopic atomic emitters in the far field (see Fig. |8|. 
For a plane-wave geometry, this would mean that the total refractive index 
(governed by the neutral and the generated plasma dispersion) of fundamen- 
tal and high-harmonic frequencies has to be identical. In typical experiments, 
HHG is carried out near the focus (waist) of a laser or in a low-diameter 
guided (wave guide, hollow fiber) geometry to realize high intensities. In this 
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case, an additional geometric wavevector component, particularly important 
for the long-wavelength fundamental beam, has to be taken into account in 
phase matching. For a focusing beam, this arises from the gradient of the 



Gouy Phase (Gouy, 1890 Siegman, 1986) (for a Gaussian beam), while for 



the guided geometry this contribution comes from the dispersion equation of 



waveguide modes (Marcatili and Schmelzer, 1964) (physically from the fact 



that the wavevector has a perpendicular component to create an effective ra- 
dial node at the waveguide wall). Phase matching in focusing beams allows 
the analysis and selection of individual microscopic quantum paths contribut- 



ing to the high harmonic spectrum (Balcou et al. , 1999). In waveguides, fully 



phase-matched HHG has been demonstrated in 1998 (Rundquist et al., 1998 



Schniirer et al. , 1998). If conventional phase matching is no longer possible, a 
spatial modulation with periodicity AL of the target medium or its environ- 
ment, resulting in an additional effective wavevector Ak = 27r/AL to balance 
the momentum conservation, can be employed to achieve quasi-phase match- 
ing (iPaul et all 120031 iBahabad et all 120101 iPopmintchev et all 120101). Also 



without a hollow-fiber guiding structure, self-guiding alone was shown to re- 



sult in efficient phase- matched HHG (|Tamaki et al. , 1999) with the potential 



for attosecond pulsed emission (Steingrube et al. 2011). The high degree of 



spatial coherence of the beam generated by a high-harmonic source has been 
demonstrated experimentally by diffraction from spatially separated pinholes 
and other objects (Bartels et al. , 2002). Below, in a), we will address the spec- 



tral coherence measurements, leading to attosecond-pulsed characteristics, in 
more detail. 



Further major breakthroughs and disruptive technological advances on the 
HHG source side included the reconstruction of the relative spectral phase 



of the high-order harmonics to prove their attosecond structure (Paul et al 



2001), the mastery and control of the carrier-envelope phase (CEP) of the 



driving laser pulses for isolated attosecond pulse generation (Baltuska et al. 



2003), the spatial and temporal control of phase matching (Pfeifer et al. 



2005bp 


i Walter et al. 


2006 


i), the gating techniques (" 


ientschel et al., 2001 


Corkum et al. 19941 


Sola et al. , 2006 Sansone et al 




2006[ 


Pfeifer et al. 


2006b|a 


i; Chang 2007 


Mashiko et al. 


2008 


Feng et al. 


2009 


|Pfeifer et al. 


2007a 


Jullien et al. 2 


008 , 


\bel et al. 


2009 


Ferrari et al. , 2010)), the com- 



prehensive control over the harmonic spectral shape towards attosecond pulse 



shaping ( 


Pfeifer et al. 




>005b 


), the creation of high-energy pulses of HHG 


light ( 


Hergott et al. 


2002 


Midorikawa et al. 


2008 


Skantzakis et al. 


2009 


Tza- 


Has et al. 


2011 


), the use of the shortest possible pulses of laser light, including 



the application of light-field synthesis methods to essentially reach the half- 



optical-cycle limit (Wirth et al. 2011), and the long-wavelength drivers to en- 



hance the high-harmonic cutoff and reduce the attosecond pulse chirp (Doumy 



et al. 2009 Popmintchev et al. 2009; Shiner et al. 2011 ). The culmination of 



this progress nowadays allows the generation of sub-100-as pulses, both iso- 



lated (Goulielmakis et al., 2008) as well as in pulse train emission mode (Ko 
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et al. , 2010), and also reaching photon energies in excess of 1 keV (Seres et al. 



2005). Advanced methods to control and selectively generate double or triple 
pulses — thus at the interface between isolated pulses and pulse trains 



arc 



currently being explored both theoretically (Raith et al 



perimentally (Pfeifer et al. 2007a Mansten et al. 2009), for applications in 



2011) as well as ex- 



attosecond interferometry. A few of these important milestones shall now be 
described in somewhat more detail. 



a) Demonstration of HHG spectral coherence and attosecond pulsed nature 



Some time after the first observation of HHG (McPherson et al. 1987; Ferray 



et al. 1988), the question arose as to whether the individual separated har- 



monic " comb" lines were coherently locked in a fixed relative phase relation. If 
their relative phases would be small on a scale of tt, the consequence would be 
an attosecond pulse train. The race was on to test this spectral phase relation. 
As attosecond pulses typically come with low pulse energies, nonlinear-optical 
autocorrelation commonly applied for femtosecond pulses was not immediately 



possible, and still is not a viable option for all photon-energy ranges (Tzallas 



et al. , 2003 Sekikawa et al. , 2004). Instead, the solution was a temporal cross 



correlation of the HHG light with a moderately intense and coherently locked 
copy of the 800 nm HHG driver pulse in a gas medium while observing photo- 



electron emission (Paul et al. , 2001). In this process, the odd harmonic photons 



/p. The ad- 



with energy E 2n +i ionized an electron with excess energy E 2n+ \ 
ditional presence of the 800 nm laser photons resulted in additional photoelec- 
tron energies at E 2n+ i — I p ±u^ by absorption or emission of 800-nm photons, 
the so-called sidebands at even-integer harmonic orders in between the peaks 
corresponding to the original odd harmonics. Importantly, interference of the 
ionization pathways proceeding via two different harmonics E 2n +i and E 2n -i 
resulting in the same final electron energies E 2n+ i — I p — uj-l = E 2n -i —I p + u^, 
allowed the detection of the relative phase between these harmonics from the 
constructive or destructive interference for the sideband peak amplitudes. This 
method has been termed RABBITT (reconstruction of attosecond beating by 
interference of two-photon transitions) and is currently used, along with its 



complementary streak-field approach (Hentschel et al. , 2001 Itatani et al. 



2002), for a variety of scientific applications (see Section |5j), including the 



measurement of photoionization dynamics. 
b) Few- cycle and CEP control technology 

After the successful demonstration of attosecond pulse trains, the next goal 
was to isolate one single attosecond pulse per driver pulse, in order to pro- 
vide a temporally singular probe for quantum-dynamics measurements. The 
group around Ferenc Krausz, being pioneers in the creation of the shortest 
optical laser pulses, were the first to measure, again by using photoelectron 



spectroscopy, an isolated attosecond pulse (Hentschel et al. , 2001). It was soon 



realized that the thus-far uncontrollable CEP of the driving laser pulse played 
a paramount role in generating single vs. double attosecond pulses per driver 
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pulse (depending on whether the driver field within the envelope exhibited 
a cos-like or sin-line shape). From a crucial cooperation with a pioneer in 
optical precision spectroscopy and creator of optical frequency comb technol- 



ogy, Theodor Hansen (physics Nobel prize in 2005 (Hansen, 2006)), the group 



was soon able to stabilize the CEP and thus the electric field within the enve- 



lope of the driver and selectively generate isolated attosecond pulses (Baltuska 



et al. 2003). Applying the same methodology with further improvements of 



femtosecond laser and optics technology, this method currently allows the pro- 



duction of the shortest isolated attosecond pulses of 80 as ( Goulielmakis et al. 



2008). 



c) HHG gating techniques 

As isolated attosecond pulses are desirable for many applications, early ideas 



focused on the isolation of attosecond pulses ( Corkum et al. , 1994 ) by means 



of employing the strong ellipticity dependence of the HHG process (Budil 



et al. 



1993|). This dependence is immediately understood from the recollision 
1993), where only linearly polarized light can lead to the 



model (Corkum 



strict classical return of the electron to the atom, while for elliptically po- 
larized light the electron would always miss its point-like origin. In quantum 
reality, it is the derealization of the electron wavefunction and the spatial 



extent of the atomic or molecular target system (Flettner et al. , 2002, 2003) 



that create a non-zero harmonic intensity response, which however remains a 
rapidly decaying function of ellipticity and typically reaches 10% of its linear- 
polarization yield for ellipticities between 0.1 and 0.2. A laser pulse exhibiting 



a temporally varying ellipticity was thus suggested (Corkum et al. 1994) and 



experimentally proven (Sola et al. 2006 Sansone et al. 2006) to result in 



isolated attosecond pulses. A major benefit of this method is the fact that the 
entire spectral bandwidth of an attosecond HHG pulse can be used for exper- 
iments, rather than only a spectral portion around its cutoff as is the case for 



linear-polarization methods (Hentschel et al. , 2001). These latter methods are 



sometimes referred to as intensity gating based on the fact that attosecond 
pulse isolation depends on the rapid variation of intensity and thus HHG cutoff 
photon energy from one half-cycle to the next within the driving laser pulse. 
The attosecond pulse has then to be selected by using a high-pass photon 
energy filter such as a metal foil before the experiment. Other gating tech- 
niques were invented that rely on two-color driving methods, i.e. two-color 

2006bp 



gating (Pfeifer et al. 



Merdji et al. 2007), also in combination with 



a time-dependent polarization field, resulting in the so-called double-optical 



gating (DOG) method (Chang, 2007 Mashiko et al. 2008), including a gen- 



eralized (GDOG) (Feng et al. , 2009) way to use pulses deep in the multi-cycle 



regime for generating isolated pulses. Yet another scheme suggests the gating 
of isolated attosecond pulses on the leading edge of the driver pulse by termi- 
nation of HHG in a rapidly ionizing medium, either by the collective-medium 



response through phase matching control (Pfeifer et al. , 2007a Jullien et al. 



2008 Abel et al. , 2009 ) or by controlling the single-atom response by com- 
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pletely depleting the ground state of the atomic target system (Ferrari et al. 
2010b. 



3 Hard x-ray HHG and zeptosecond pulses 



The extension of HHG sources toward higher photon energies and shorter 
pulse durations is cr ucial to enable new applications such as, e.g., time-resolved 
diffraction imaging (Neutze et al. , 2000) with sub- Angstrom resolution, excita- 



tion of tightly bound core electrons or even the time-resolved study of nuclear 



excitations (Burvenich et al. 2006 Weidenmuller , 2011). The advancements 



in the time and energy domain are inter-connected. A larger harmonic band- 
width is required to generate shorter light pulses. Presently, the keV regime 
of HHG is approached allowing to break the one attosecond barrier toward 
pulses of zeptoseconds duration as visible from the time-bandwidth product 
[Eq. d9b 



The spectral cutoff energy of HHG [see Eq. (pj)] can principally be increased 
by employment of higher laser intensities. This way, photon energies of few 



keVs have been generated (Seres et al.l 2005, 2006) with small rates just above 



the detection threshold [see Fig. |3j (b)] because this straightforward increase 
of the laser intensity is connected with several difficulties. First, rising the 
laser intensity above resulting harmonic cutoff energies of few hundred elec- 
tronvolts leads to a large free-electron background causing phase mismatch 
and tiny emission yields on a macroscopic level. For this reason, in a different 



approach ( 


Popmintchev et al. , 


2008 


Doumy et al. 


2009 


Popmintchev et al. 


2009 


Chen et al. 


2010 


Ghimire et al. 


2011 


), higher cutoff energies are ob- 



tained without significant change of the free-electron background by increase 
of the driving laser wavelength for fixed laser intensity. This way, increasing 
the ponderomotive potential without increasing intensity, one can potentially 
attain photon energies corresponding to the nonrelativistic HHG limit. Sec- 
ond, relativistic effects can lead to a dramatic suppression of the single atom 



HHG yield (Di Piazza et al. 2012). While the nonrelativistic electron mo- 
tion is mainly along the polarization direction of the laser field [see Fig [9^,] , 
the Lorentz force becomes noticeable for electron velocities approaching the 
speed of light. Then the electron is transfered along the propagation direction 
of the laser pulse hindering recombination and frustrating HHG, see. Fig |9Jd. 
The relativistic drift becomes significant if the drift distance exceeds the wave 
packet dimensions after spreading. This happens at laser intensities of about 
10 17 W/cm 2 ( Palaniyappan et al. 2006) at 800 nm wavelength and marks the 



limits in terms of photon energy for the nonrelativistic HHG at about 10 keV. 



For short pulse generation, besides a large bandwidth, synchronization of the 
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Fig. 9. (Color online) Two different classical electron trajectories for a linearly po- 
larized laser field in the (a) nonrelativistic and (b) highly relativistic regime. In the 
relativistic case the electron does not return to the atomic core when the laser field 
changes its sign. HHG is unlikely what is depicted by the dashed line. 

emission of the different harmonic components in the pulse is required. How- 



ever, the HHG pulses have an intrinsic chirp, the so-called attochirp (Mairesse 



et al. , 2003 Doumy et al. , 2009 ) and, thus, are much longer than their band- 



width limit. To compress the emitted pulse down to its fundamental limit, 



dispersive elements of either chirped multilayer x-ray mirrors (Morlens et al 



2005), thin metallic films (Kim et al. , 2004; Lopez-Martens et al. , 2005), grat 



ing compressors (Poletto et al. , 2008) or thick gaseous media (Ko et al. , 2010 



Kim et al. 2007) are employed. However, with increasing bandwidths, it will 



be difficult to find media with suitable dispersion in the future. 



In the following a selection of techniques for enhancing the cutoff energy or 
decreasing the pulse duration of the harmonic light are reviewed. 



3.1 HHG with long wavelength drivers 



As indicated in the previous part, the quest of reaching keV photon energies 
and zeptosecond pulses is ultimately linked with solving the phase-matching 
problem. For coherent growth of the harmonic signal, the phase velocities of 
the harmonic and laser field are required to be equal. Under this condition, the 
light emitted from different atoms adds up constructively as sketched in Fig. [8] 
and explained in Sec. |2.4| The condition of equal phase velocities can also be 
defined as a vanishing phase- mismatch vector Ak (Gaarde et al. , 2008). The 



phase-mismatch vector is (Durfee et al. 1999 Popmintchev et al. 2009) 



Ak « qk g 



- qp(l - riy^-iAb + n 2 ) + qpqN^rcXi, + a^- 
Al dx 



(13) 



where q is the harmonic order, p the gas pressure, Al the laser wavelength, 
AS the difference of the index of refraction per pressure between the laser 
wavelength and the harmonic wavelength, n 2 = n 2 Ih the nonlinear index of 
refraction per pressure, iV a the number density of atoms per pressure, and r e 
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the classical electron radius. qk geo is a geometric term due to waveguide or 
near-focus propagation. The second and third term cover the refractive index 
difference between both waves caused by atomic and plasma dispersion, re- 
spectively. The last term covers the variation of the atomic dipole emission 



phase (Lewenstein et al. 1995) with the laser intensity introduced in Eq. (12). 
Moreover, non-adiabatic effects due to the laser pulse deformation, which are 
important in the case of ultra-short driving pulses especially close to the sat- 
uration regime, are excluded, see, e.g., Tempea et al. (2000); Geissler et al. 



(2000). For a detailed review on the phase-matching aspects of HHG see, 



e.g., Gaarde et al. (2008). 



In experiments, conditions have to be found to achieve vanishing Ak [Eq. ( 13 )] . 



We concentrate on the discussion of waveguide propagation where a ^ is negli- 
gible. It is beneficial that the two pressure-dependent terms of the atomic and 



plasma dispersion in Eq. (13) have opposite signs and, thus, can cancel each 
other to some extent for low ionization degrees r\. Hence, phase-matching is 
commonly achieved by selecting the pressure p such that Ak vanishes exploit- 
ing that the pressure independent geometric term qk geo balances the difference 
between the two dispersion terms that depends linear on p. At a fixed IR laser 
wavelength but increasing laser intensity, balancing is possible to some extent 
to attain HHG energies up to about 100 eV depending on the atomic species. 
For HHG cutoff energies above this limit, the laser intensity is sufficiently 
high to ionize a considerable part r\ of the atoms leading to a serious misbal- 
ance between the atomic and plasma dispersion. Above this critical level of 
ionization rj CT , it is not possible to counteract the misbalance by an increase 
of the gas pressure any more. An obvious approach to circumvent this prob- 
lem is to keep the laser intensity on a low level but to increase the driving 
wavelength. This way, the quadratic scaling of the single-atom cutoff energy 
[Eq. (J7])] with the driving wavelength is exploited while the ionization level 
r\ remains constant because the ionization rate in the tunneling regime de- 
pends mainly on the laser intensity rather than its wavelength. Long time, the 
approach was not pursued because the single-atom yield has an unfavorable 
scaling with the driving wavelength of 



A 



-5.5±0.5 



(Corkum 



1993 



Gordon 



and Kartner 2005a 


Tate et al. 


, 2007 


( 


Colosimo et al. 


2008 


Frolov et al. 


2008[ 


Doumy et al. |2009 ; Schiessl et al. 


2007 


Perez- Hernandez et al. 


2009 



Shiner et al. 



2009). Despite this, the generation of bright soft x-ray harmon- 



ics with mid-infrared drivers has been demonstrated recently (Popmintchev 



et al. 2008, 2009 Chen et al. 2010) exploiting the fact that near the critical 
level of ionization, high gas pressures are required to achieve phase-matching. 
For high gas pressures, the increased number of emitters compensates the de- 
crease of the single-atom response yielding efficient overall HHG. Gas-filled 



capillaries are especially suitable in this regard (Durfee et al. , 1999). This way 



efficient generation of HHG with several hundreds of electronvolts bandwidth 



has been demonstrated ( Chen et al. , 2010 ) being sufficient to support pulses 



of only lias duration. From this point of view, the concept of increasing the 
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driver wavelength seems highly promising at the moment. 



In the case that the phase-matching condition Ak « is not feasible, so-called 
quasi-phase- matching (QPM) schemes are invoked. Their principle relies on 
symmetry breaking between the positive and negative contributions responsi- 
ble for the destructive interference. In this regard a weak counterpropagating 



IR field (Peatross et al. , 1997; Cohen et al. 2007), weak static fields (Serrat 



and Biegert, 2010), modulated wave guides (Christov et al. 



2003 Gibson et al. 2003) or multiple gas jets (Willner et al. 



2000 



Paul et al. 



2011) can be em- 



ployed. However, these schemes require additional efforts in the experimental 
implementation. 



3.2 Relativistic regime of HHG 



Although the approach in Sec. 3.1 to increase the driver wavelength raises 
the perspective to construct bright HHG sources above one keV, limits will 
appear in the multi-keV regime for several reasons. First, the velocity of the 
electron will become comparable to the speed of light and, thus, the electron 
will undergo a drift motion caused by the Lorentz force preventing recollision 
along with recombination [see also the discussion in the beginning of Sec. [3] . 
Second, the increase of the driver wavelength is accompanied by a longer 
electron excursion which can exceed the mean distance between the atoms at 
the required high gas pressures. The resulting collisions of the electron with 



other atoms are a source of incoherence and prevent efficient HHG (Strelkov 



et al. 2005). Facing these difficulties, strategies are required to overcome the 
relativistic drift and to solve the phase-matching problem in a gas of multiply- 
charged ions rather than using long-wavelength drivers on a low ionization 
level. 



While in the weakly relativistic regime the drift problem is not severe and 



harmonic emission can be observed (Latinne et al. , 1994 Hu and Keitel, 2001 



Keitel and Hu, 2002), at higher intensities special methods have to be applied 



to counteract the relativistic drift. To suppress the drift, highly charged ions 



moving relativistically against the laser propagation direction (Mocken and 



Keitel 


, 2004 


et al. , 


2006) 


et al. 


2006) 



Chirila et al. , 2004), antisymmetric molecular orbitals (Fischer 



2006) or a gas of positronium atoms (Henrich et al. , 2004 Hatsagortsyan 



2006 ) can be used. Different combinations of laser fields have also been 



proposed for this purpose such as a tightly focused laser beam (Lin et al 



2006), two counterpropagating laser beams with linear polarization (Keitel 



et al. 1993 Kylstra et al. 2000 Taranukhin, 2000 Verschl and Keitel, 2007b) 



or with equal-handed circular polarization (Milosevic, N. et al. , 2004 ; Verschl 



and Keitel, 2007b). In the latter field configuration, the relativistic drift is 



eliminated. However, in this scheme the phase-matching is particularly prob- 
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lematic to realize (Liu et al. 2009). In the weakly relativistic regime, the 



Lorentz force can also be compensated by a second weak laser beam being 



polarized in the strong beam propagation direction (Chirila et al. , 2002). Two 



consecutive laser pulses or the laser field assisted by a strong magnetic field 
have been proposed as well (Verschl and Keitel, 2007a cl) but this requires 



large magnetic fields and dilute samples. In Klaiber et al. 



2006, 2007a) it has 



been shown that the relativistic drift can be significantly reduced by means 
of a special tailoring of the driving laser pulses in the form of attosecond 
pulse trains (APTs). These results concern the relativistic drift problem only. 
However, in order to generate relativistic harmonics, both the drift and the 
phase-matching problems have to be tackled simultaneously. Below we re- 



port on two works (Kohler et al. 2011b Kohler and Hatsagortsyan, 2012) 



approaching those jointly. 



In Kohler et al. (2011b), HHG emission from a gaseous medium driven by 



two counterpropagating APTs has been investigated. The field configuration 



is capable of circumventing the impact of the relativistic drift (Klaiber et al 



2007b Hatsagortsyan et al. 2008) as seen from the schematic of the electron 



dynamics in Fig. 10 (a). The electron is liberated by laser pulse 1 which reaches 



the atom first. Subsequently, it is driven by this pulse in the continuum and 
undergoes the relativistic drift. This part of the trajectory is indicated by the 
light blue coloring [in Fig. [To] (a)]. Thereafter, the electron propagates freely 
(gray dashed) in the continuum, before, a moment later, the second pulse 
reaches the electron, reverts the drift and imposes recollision (dark blue). This 
setup has a distinguished property which facilitates the realization of phase 
matching despite the significant free-electron dispersion. In fact, for different 
atoms situated along the propagation direction, the time delays between the 
first pulse and the second pulse are different. This is reflected in different 
intrinsic harmonic emission phases. The phase can be tuned by variation of 
the laser field intensity to compensate the phase mismatch caused by the free 
electron background. The period of the spatial modulations is given by the 
APT period. 



In Klaiber et al. (2008); Hatsagortsyan et al. (2008); Kohler and Hatsagortsyan 



(2012) relativistic HHG is achieved employing a setup where the laser field is 



assisted by an x-ray field. The usefulness and applicability of x-ray assistance 
have been proved in various theoretical investigations and experiments for 
different purposes (see Sec. 3.3). In the present case, the assisting x-ray field 



is required for the ionization via one x-ray photon absorption. Therefore, the 
x-ray frequency has to exceed the ionization energy conveying the electron 
a large enough initial momentum opposite to the laser propagation direction 
which is able to compensate for the relativistic drift. The setup is displayed 



in Fig. 10 (b). A collinear alignment of the laser and x-ray field propagation 



directions is shown to be advantageous in terms of phase-matching (Kohler 



and Hatsagortsyan, 2012). An additional weak counterpropagating IR field 
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Fig. 10. (Color online) (a) Single-atom perspective of the HHG setup with coun- 
terpropagating APTs. The classical trajectory of a rescattered electron of a single 
atom in the gas target. After ionization by pulse 1, the ejected electron is driven 
in the same pulse (light blue, light gray), propagates freely after the pulse has left 
(gray dashed) and is driven back to the ion by the second laser pulse (dark blue, 



dark gray). Figure reprinted with permission from Kohler et al. (2011b). Copyright 
2011 by the Institute of Physics, (b) Geometry of the HHG process for a collinear 
alignment of the x-ray and laser fields. A weak IR field is added that accomplishes 
phase- matching. 

arranges a quasi-phase- matching scheme. 

Both of the former relativistic HHG setups provide a tiny macroscopic HHG 
yield of 10 photons per shot for a HHG energy of about 50 keV in common. 
Several issues have been identified lowering the HHG yield when extending 



the HHG cutoff towards the multi keV regime (Kohler et al. , 2011b Kohler 



and Hatsagortsyan, 2012). To discuss them, the macroscopic yield can be 



approximately expressed by the product 



N 



2 „2 



x An x Att x At x V z p 



(14) 



where dw n /dQ is the single-atom emission rate, An is the number of harmon- 
ics within the phase-matched frequency bandwidth, At the interaction time 
that is approximately the delay between both pulses, Aft the solid angle of 
emitted harmonics, V the volume of coherently emitting atoms (perfect phase- 
matching is assumed in this volume) and p the atomic density. The terms in 



Eq. ( 14 ) can be evaluated for the two relativistic setups as well as for an exper- 



iment with standard HHG in the nonrelativistic regime, see e.g., Goulielmakis 
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et al. (2008). Comparing the factors in Eq. (14) between the different setups 



one notes that in the relativistic regime suppression arises from the single-atom 
yield, the smaller solid emission angle and the lower gas density. The reasons 
for the lower single-atom yield are manifold: The higher momentum of the 
recolliding electron favors scattering instead of recombination. Furthermore, 
the larger energy range of the recolliding wave packet leads to a smaller HHG 
yield per harmonic. Then compared to a non-relativistic sinusoidal field only, 
the ionization and continuum dynamics of the wave packet are changed and a 
smaller fraction actually revisits the core area. In analogy to the interference 
pattern of an aperture, the solid emission angle decreases quadratically with 
the wavelength. Additionally, phase-matching could only be accomplished for 
either low density (10 16 cm -3 in Kohler and Hatsagortsyan (2012)) or small 



volume (length ~ 12.5 /xm in Kohler et al. (2011b)). Note that the formerly 
discussed reasons of surpression are not caused by relativistic effects rather 
than by the high HHG energy in general. 



On the bottom line, the former results show that more efforts are required to 
improve the efficiency of HHG in the multi keV regime. Below we will discuss 
another way to increase the HHG energy without increasing the laser intensity 
and suffering from drift problems. 



3.3 XUV-assisted HHG 



The usefulness of XUV light assisting a strong laser field has been demon- 
strated for various purposes. It has been used to enhance HHG by many orders 



of magnitude compared with the case of the laser field alone (Ishikawa, 2003 



Takahashi et al. , 2007 ) . When the XUV field has the form of an attosecond 
pulse train a single quantum path can be selected to contribute to HHG and in 
this way allows to manipulate the time-frequency properties of the harmonics 



as well as to enhance a selected bandwidth of harmonics ( 


Schafer et al. , 


2004 


Gaarde et al. 


2005 




?igueira de Morisson Faria and Salieres 


2007 


), and to 


extend the spectrum 


{ Fleischer and Moiseyev 


2008 


)• 





Tuning an intense XUV field from an FEL to a resonance between a core and 
valence state can lead to the emergence of a second plateau that is shifted 
to higher energies by the former resonance energy with respect to the first 



plateau (Buth et al. , 2011). The presence of originally two bound electrons is 



thereby crucial for the effect. Plateau extension in the presence of two electrons 



of the secondary plateau. In Buth et al. 



has already been noticed in Koval et al. (2007), however, with a low probability 



(2011), the intensity of the secondary 



plateau is tunable via the FEL intensity. A schematic of the proposed scheme 



is shown in Fig. 11 The atoms are irradiated by both an intense optical laser 
field and the resonant x-ray field from an FEL. As soon as the valence electron 
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(a) 



(b) 




Fig. 11. (Color online) Schematic of the HHG scenario as a three-step process: (a) 
the valence electron is tunnel ionized; (b) the additional high-frequency light excites 
the core electron; (c) the continuum electron recombines with the core hole. Figure 



reprinted with permission from Buth et al. (2011). Copyright 2011 by the Optical 
Society of America. 

is tunnel ionized by the optical laser field, the core electron can be excited to 
the valence vacancy by the x rays. Then the continuum electron, returning 
after a typical time of 1 fs, can recombine with a core hole rather than with 
the valence hole from that it was previously tunnel ionized and thus emit a 
much higher energy. 



In Buth et al. (2011 2012) an analytical formalism is developed to cope with 
the two-electron two-color problem. Losses due to tunnel ionization and direct 
x-ray ionization are included via phenomenological decay constants in con- 
junction with Auger decay of the intermediate hole. The theory is applied to 
the 3d — >• 4p resonance in a krypton cation as well as to the Is — > 2p reso- 
nance in a neon cation. The results for a resonant sinusoidal x-ray field for 



two different intensities are shown in Fig. 12 The chosen optical laser field 
intensity is 3 x 10 14 W/cm 2 for krypton and 5 x 10 14 W/cm 2 for neon both at 
800 nm wavelength. 
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Fig. 12. (Color online) HHG photon numbers of the n-th harmonic for different x-ray 
intensities from (a), (c) krypton and (b), (d) neon. The solid black line stands for the 
recombination to the valence state whereas the more energetic plateau arising from 
core hole recombination is the red dashed line. The XUV and laser pulse durations 



are three optical laser cycles in all cases. Figure adapted from Kohler et al. (2012). 



The most striking feature in the obtained spectra is the appearance of a 
second plateau. It is up-shifted in energy with respect to the first plateau 
by the energy difference between the two involved core and valence states. 
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The two plateaus have comparable harmonic yields for x-ray intensities above 
10 16 W/cm 2 . Quite importantly, the losses due to x-ray ionization do not lead 
to a significant drop of the HHG rate. The second plateau bears signatures 
of the core state and may offer a route for ultrafast time-dependent chemical 
imaging of inner shells (Itatani et al. 2004 Morishita et al. 2008). Moreover, 
by exploiting the upshift in energy, attosecond x-ray pulses come into reach. 



3.4 Exotic light sources 



A route to shift the attainable cutoff energy to the multi-keV regime with- 
out suffering from the relativistic drift is the use of muons rather than elec- 
trons bound to the atomic system. According to the higher mass of the muon 
« 207 compared to the electron mass m e = 1, its kinetic energy is in- 
creased by a factor of the reduced mass m r = m^nin/ (m M + m n ) for the same 
velocity where m n is the nuclear mass. Thus, the attained energy of a muon is 
much higher for nonrelativistic velocities than for electronic systems. This can 
be expressed more quantitatively by the linear scaling of the ponderomotive 
energy U p = at the non-relativistic limit defined by the field strength 

parameter £ rel which has a weak dependence on m r ( Hatsagortsyan et al 



2008). Note that higher laser intensities are required to accelerate the muon 
to the non-relativistic limit as can be seen from £ = — - — . The high recolli- 



sion energies could allow to trigger nuclear reactions ( |Chelkowski et al~ 2004), 
to probe the nuclear structure due to the smaller Bohr radius of the bound 



muon (Shahbaz et al., 2007) and to create zeptosecond pulses (Xiang et al. 



2010). A precaution for this kind of experiments is certainly the realization 



of a target density being sufficiently high to generate a measurable emission 
yield. 



Based on a different approach, Ipp et al. (2009) proposed a method devised 



to generate yoctosecond pulses. The authors show that quark gluon plasmas 
generated by high-energy heavy ion collision can be a source of flashes of GeV 
7-rays having such a short duration. Thereby, the fast expansion of the quark 
gluon plasma after collision is responsible for the short time scale. Under cer- 
tain conditions a momentum anisotropy of the expanding plasma may form 
at intermediate times which leads to a preferential photon emission direction 
perpendicular to the collision axis during that short time. Therefore, a double 
pulse structure is emitted then in the other directions which has a princi- 
pally variable pulse delay. The scheme could be employed for pump-probe 
measurements on the yoctosecond time scale. 
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4 HHG in shaped driving pulses 



Femtosecond pulse shaping of the driving laser field (Weiner et al. , 1990; Zeek 



et al. 1999 Kornaszewski et al. 2008) is one of the most direct ways for 



controlling HHG. It allows for the precise adjustment of the ionization rates 
and the control of the electron continuum dynamics and, thus, for almost 
arbitrarily engineering the recolliding electronic wave packet with impact on 
the spectral and temporal properties of the emitted harmonic light. 



4-1 HHG yield and cutoff enhancement 



First experimental results using shaped driver laser pulses have demonstrated 
the enhancement of HHG, including the selective enhancement of individual 



harmonic orders (Bartels et al. , 2000). Using similar methods, control over the 



spectral positions of the produced harmonic spectral comb lines has also been 



demonstrated (Reitze et al. 2004). 



On the theoretical side, Chipperfield et al. (2009) have recently proposed an 



optimal pulse shape for reaching the maximum cutoff under a given laser pulse 



energy and laser period. The saw-tooth-form with a DC offset (red) in Fig. 13 1 



turned out to be optimal (Radnor et al. , 2008 Chipperfield et al. , 2010) with a 



cutoff energy about 3 times higher than a sinusoidal pulse of the same energy. 
The respective classical trajectory leading to the cutoff energy is shown as 
green dashed line. The optimal waveform has its highest intensity at the end 
of the cycle because this way it nurtures the photon energy the most just 
before recollision. The electric field in the initial part of the trajectory is only 
required to drive the electron away. High intensities along with high velocities 
in this part of the trajectory are not required for reaching the maximum 
cutoff but are at the expense of the pulse energy and, thus, would reduce the 
electric field required for the final acceleration. For this reason the intensity 
is lower in the beginning of the trajectory than at the end. In Fig. 13: the 
authors demonstrate that the first four Fourier components are sufficient to 
reproduce the optimal waveform to an adequate level. However, the pulse form 
lacks in efficiency with regard to the generated harmonics because the initial 
electric field is only about half of the final electric field. Thus, depending on 
the ionization potential the ionization rate is either small or over-the-barrier 
ionization happens at the end of the pulse suppressing HHG (see Sec. 2.2). 



For this reason, the field is modified in the beginning of the period by a field 
of doubled wavelength to enhance the ionization rate. Quantum simulations 
demonstrate a similar efficiency for the generated harmonics compared to the 
case of a sinusoidal field. 
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Fig. 13. (Color online) The solid red lines represent the electric field whereas the 
green dashed line is the respective electron excursion, (a) Optimal waveform in 
terms of recollision energy, (b) Sinusoidal waveform similar to the form in (a) but 
with a longer period (1.54T). (c) Waveform of (a) composed of the first four Fourier 
components, (d) As (c), but optimized with respect to the ionization rate. Figure 



reprinted with permission from Chipperfield et al. (2009). Copyright 2009 by the 
American Physical Society. 

Pulse shaping can also be employed to suppress the relativistic drift. The 
classical trajectory of a recolliding electron can be divided into a part where 
the electron moves away from the core and another where it returns to the core. 
In both parts of the pulse, the electron undergoes a relativistic drift motion. To 
reach the final high relativistic recollision energy the electron is not required 
to move during the first part of the trajectory with relativistic velocities in the 



strong field. For that reason the pulse is shaped in Klaiber et al. (2006) such 
that the field is negligible during the first part of the trajectory and, thus, the 
drift is avoided then. 



4-2 Attosecond pulse shaping 



Full control over attosecond pulse shapes however requires the comprehen- 
sive control over the HHG spectrum in both its amplitude and phase. First 
experiments along these lines have focused on the control of the broadband 
harmonic spectral shape (Pfeifer et al. , 2005b), suggesting a collective medium 



response (phase-matching) mechanism for the observed controllability (Pfeifer 
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et al. 2005a Walter et al. 2006). Groups of harmonics as well as various 



individual harmonics could be selectively generated. Also suppression of a 
single harmonic peak in an extended spectrum was possible. The former 
spectral-amplitude shaping capability alone can result in major modifications 
of the attosecond pulse shapes that are produced. An experiment applying 
such variably shaped light fields to an SFg molecular target system in an 
approach towards coherent electron control with shaped attosecond pulses 
demonstrated selectivity in the branching ratio of dissociative photoionization 



channels (Pfeifer et al. 2007b). Recent advances in the application of pulse 



shaping to HHG have also been reviewed in Winterfeldt et al. (2008). 



The former paragraph mainly dealt with the engineering of the spectral HHG 
amplitude. Phase control of high-harmonic spectra has been studied by means 



of passive dispersive effects (Lopez-Martens et al. 2005 Strasser et al. , 2006) 



or by active laser control of molecular alignment (Boutu et al. , 2008). A multi- 



dimensional control scheme including the CEP and broadband spectral driver 
pulses in the few-cycle regime for enhanced attosecond pulse control of am- 
plitude and phase has also recently been suggested (Raith et al. , 2011). Once 



available, fully controlled and intense attosecond fields may enable applications 
such as the preparation of exotic atomic and molecular electronic states, pos- 
sibly leading to the creation of novel types of bonding which is out of reach 
of traditional (thermodynamics) methods. It may also provide access to an 
unimaginable technology in molecular electronics, where (possibly quantum) 
information is created upon, processed and computed while traveling along 
and crossing tiny molecular wires exhibiting conjugated electronic bonds and 
the correspondingly delocalized electronic wavefunctions. 



Nowadays pulses down to a duration of 63 as (Ko et al. , 2010) have been 



generated and the bandwidth to create pulses of only 11 as is available (Chen 



et al. , 2010). Although the emitted harmonic light possesses a large bandwidth 



its time structure does not automatically correspond to its Fourier limit [see 
Eq. ([9])] because of the attochirp stemming from the HHG process. To com- 
press the emitted pulse down to its fundamental limit, dispersive elements are 
employed as mentioned earlier. However, it may become demanding in the 
future to find suitable elements for multi keV bandwidths. An alternative way 
to circumvent the attochirp problem would be to modify the HHG process 
such that the light is emitted without attochirp. To achieve this goal, the con- 



tinuum dynamics can be altered by adding a second-harmonic (Zheng et al 



2009) or a subharmonic (Zou et al. , 2010) to the laser field. Significant chirp 



compensation has been proved in this case. 



A method to acquire a complete control of the chirp has been proposed 
in Kohler et al. (2011a) by means of laser pulse shaping and soft x-ray as- 



sistance using an ionic gas medium. This method allows for the formation of 
attosecond pulses with arbitrary chirp, including the possibility of attochirp- 
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n/2 
Laser phase 

Fig. 14. (Color online) Schematic of the recollision scenario: a) A half cycle of the 
tailored laser field (black). The wiggled red (gray) line is the field of the assisting 
x-ray pulse, b) Different one-dimensional classical trajectories in the field of (a) 
which start into the continuum at different times but revisit the ionic core at the 



same time. Figure reprinted with permission from Kohler et al. (2011a). Copyright 
2011 by the Optical Society of America. 



free HHG and bandwidth-limited attosecond pulses. The principle is illus- 
trated in Fig. [14] which shows the relevant trajectories in the shaped field 



providing HHG. Since the recollision time of a certain harmonic can be iden- 



tified with its group delay (Mairesse et al. , 2003; Kazamias and Balcou, 2004) 



in the emitted pulse, a simultaneous recollision of all trajectories leads to a 
bandwidth-limited attosecond pulse. The demand of simultaneous recollision 
can be fulfilled if the electron is freed by single-photon ionization when the 
x-ray frequency wx is much larger than the binding energy. In this case, the 
electron has a large initial kinetic energy directly after ionization. Let us focus 
on the two example trajectories marked by a and f3 in Fig. [Tip . Both are ion- 
ized at instants separated by a small time difference 5t\. Because a experiences 
the laser field a little earlier, a velocity difference arises between both. Note 
that the velocity difference between a and is conserved in time for a homo- 
geneous laser field. With a convenient choice of the parameters, the velocity 
difference between both acquired during 5t\ can be such that both recollide 
simultaneously. In principle, the single-photon ionization process allows for a 
variety of different ionization directions. However, only the trajectories start- 
ing exactly opposite to the polarization direction of the laser recollide and are 
thus considered. 



The laser field with a shape approaching the ideal one shown in Fig. |14fi 
can be obtained using a small number of Fourier components. The required 
parameters which are sufficient to reach the bandwidth-limited pulses below 
10 as and 1 as are indicated in Table [TJ However, the method still has several 
drawbacks. It is experimentally demanding to create a pure ionic gas and 
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N F 


/ L [W/cm 2 ] 


w x [eV] 


/ x [W/cm 2 ] 


ion 


I P [a.u] 


At [as] 


8 


10 16 


218 


3.5 x 10 14 


Li 2 + 


4.5 


8 


20 


10 17 


996 


1.4 x 10 15 


Be 3+ 


8 


0.8 



Table 1 

Np represents the number of Fourier components contained in the tailored driving 
pulse, Jl its peak intensity, wx the x-ray frequency employed for ionization, 7x its 
intensity, I p the ionization energy and At the harmonic pulse duration. 

to achieve phase-matching in a macroscopic medium due to the free-electron 
dispersion. Moreover, the ionization rate is small due to wx 3> I P , and the 
required large initial momentum and the dipole angular distribution of the 
ionization process lead to an increased spread of the ionized wave packet as 
compared to tunnel ionization. Additionally, precise shaping of intense driver 
fields with intense harmonics as well as the synchronization of the x-ray and 
IR pulses is also experimentally demanding. 



5 Experimental applications 



Continuously gaining momentum from the increased knowledge on HHG and 
the co-evolution of technological capabilities, the scientific applications of at- 
tosecond HHG light have begun and keep radiating out into a range of various 
research directions. Owing to their short-pulsed nature and high photon en- 
ergies available for probing, the old scientific dream of tracking and watching 
the dynamics of individual and multiple electrons on their quantum paths 
throughout the valence and core-shells of atoms and molecules was rekindled. 
Using HHG-based sources, the scientific community has already begun to ob- 



serve the first valence-shell wavefunction dynamics in atoms (Goulielmakis 



et al. 2010) and molecules (Smirnova et al. 2009), and will soon explore 



wavepackets of two (Argenti and Lindroth, 2010) or multiple electrons such 
as in plasmon excitations (Pfeifer et al. 2008 Siifimann and Kling, 2011) in 
motion. While the femtosecond revolution already opened up the now well- 
established and fruitful field of femtochemistry (Ahmed Zewail, Nobel prize 
1999 (Zewail, 2000)) making intramolecular motion of atoms accessible to 
experimentalists, attosecond pulses now promise to be key to monitoring in- 



dividual bonds in molecules being made or broken ( Nugent- Glandorf et al. 



2001 Wernet et al. 2009; Worner et al, 2010), directly clock the instant of 



photoionization (Cavalieri et al. , 2007 Eckle et al. , 2008 Schultze et al. 2010) 



or to follow electronic excitations throughout small molecules (Remacle and 



Levine, 2006). The experimental capabilities in terms of temporal resolution 



have come far enough to ask truly fundamental physical questions: Are the 
physical concepts we use for defining "times" (such as ionization or tunneling 
time) valid? Are some questions about the "timing" of electronic processes 
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even ill posed and rather addressed by the definition of phases? How does 
the Coulomb potential surrounding atomic and molecular ions, modify the 
interpretation of measured timings? And to what extent does the presence of 
a moderately strong laser field, which typically is involved in measurements, 
modify the timing of some such fundamental processes? Here, we will outline 
a few selected, and by no means exhaustive, frontiers of current experimen- 
tal approaches invented to measure attosecond electronic quantum dynamics. 
What is common to all these methods is the now old-fashioned concept of 
pump-probe spectroscopy: a temporally compact trigger event (pump), which 
can be an attosecond pulse or an electric-field extremum within an intense 
optical laser cycle and a similarly generated probe event, which can either be 
controlled in time or the temporal position of which can be measured. 



a) Photoelectron spectroscopy methods for time-resolving ionization dynamics 
The first and still most widespread experimental method for time-resolving 
attosecond processes is the detection of photoelectrons. Owing to its typical 
photon energy in the vacuum ultraviolet (VUV) up to the soft-x-ray region, 
HHG-light interaction with neutral matter almost exclusively results in ion- 
ization of the system under study. The first experiments (see Section 2.4) 
on the characterization of attosecond pulses proceeded through the observa- 
tion of energy-resolved photoelectrons as a function of delay time between the 
attosecond and a coherently locked femtosecond pulses, both in the form of 
attosecond pulse trains (RABBITT (Paul et al. 2001)) or isolated attosec- 
ond pulses (streaking (Hentschel et al. 2001 Itatani et al. 2002)). While 
for the RABBITT technique, the so-called atomic phase arising from single- 
photon ionization (a system- and initial-state specific photon-energy uj depen- 
dent phase ip(uj) of the bound-free dipole matrix element) was always included 
to interpret the results, the streaking technique initially assumed an "instan- 
taneous" response (At = d< ^^ ~ 0) of the photoionization process, i.e. the 
working principle that the photoionization rate follows the exact intensity 
shape of the attosecond pulse. Also due to this assumption, experimental re- 
sults on the ionization of solids (Cavalieri et al. , 2007) came as a surprise, 
where it was noticed, by streak-field cross-correlation of a 91 eV attosecond 
photoionizing pulse with a moderately intense near-visible pulse that elec- 
trons seem to ionize "sooner" or "later" depending on their initial state. The 
photoionization time difference between the two initial states being either in 
the valence band or a more localized 4d inner-valence state was measured 
to be 110±70 as, and tentatively interpreted as being caused by the differ- 
ent band structure and different transport processes of the electrons on their 
way to the surface, where they are further accelerated or decelerated by the 
near-visible laser field for streak-field detection. Other experiments followed 
on the photoionization of gas-phase atoms, where similar temporal shifts of 
photoionization times were discovered. For Ne ionization by an attosecond 
pulse, a 20 as "delayed emission" of the 2p vs. the 2s states was found by 
the same experimental method and tentatively discussed in the context of 
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electron correlation (Schultze et al. 2010). As the strong-field approximation 



was used in reconstructing the temporal shifts, and also as the moderately 



intense laser pulse modifies the dynamics (Kliinder et al. , 2011; Ivanov and 



Smirnova, 2011), these and the previous solid-state results are still subject of 



active and controversial discussions 


' Lemell et al. 


2009; 


Zhang and Thumm 


2009 


Baggesen and Madsen, 


2010 : 


Klreifets and Ivanov, 


2010| 


Nagele et al. 


2011 


. It is certain, however, that this experimental approach has reachec 



a level of sophistication that tiny temporal or phase shifts can be extracted, 
where it is now up to advanced theory to interpret the results of such measure- 
ments by inclusion of the Coulomb potential, the probing laser field and the 



already well-known Eisenbud-Wigner- Smith time delay (Ivanov and Smirnova 



2011 Wigner, 1955 Smith, 1960), to extract additional knowledge on atomic 



or molecular (Haessler et al. , 2009 Caillat et al. 2011) systems beyond the 



single-active electron picture. 

b ) HHG recollision spectroscopy for measuring electronic wavepackets in molecules 
Instead of measuring photoelectrons or photoions produced in the interaction 
of HHG attosecond pulses with matter, the HHG process itself can yield dy- 
namical information about electron motion within quantum systems. With 
traditional HHG being due to interference of the laser-accelerated recolliding 



continuum states with the bound initial state of the electron (see Section 2.2), 
the spectrum of the HHG light contains information, even dynamical infor- 
mation, on the shape of the initial molecular orbital ionized by the strong 
laser field. By analyzing HHG from impulsively (field-free) aligned molecules 
and by comparison of the alignment-angle-dependent HHG spectra to a suit- 
able atomic reference target, the highest-occupied-molecular-orbital (HOMO) 
shape of nitrogen was reconstructed in a tomographic way ( Itatani et al.| 
2004). Already in this first work, the potential of this method with respect to 



observing electronic wavepacket dynamics was pointed out. After first signs 
of the contribution of several molecular orbitals (e.g. HOMO and HOMO- 



1) to the HHG process (McFarland et al. 2008), coherent interference be- 
tween these channels was measured in amplitude and phase and analyzed 
in dependence of alignment angle, HHG photon energy and driving laser in- 
tensity. From this analysis, the wavefunction shape of the HOMO/HOMO-2 
superposition at the time of tunnel-ionization in the HHG process was iden- 



tified (Smirnova et al. 2009). In these experiments, one makes use of the 
fact that time-resolved information is encoded in the harmonic spectrum, as a 
certain recolliding (probe) electron energy (and thus emitted photon energy) 
corresponds to a certain time-delay after the tunnel-ionization (pump) by an 
electric field maximum. Additional measurement of the ellipticity of the HHG 
light as a function of alignment angle and photon energy recently revealed 
information on the dynamical hole motion after strong-field ionization of ni- 
trogen molecules, including sub-cycle inter-orbital population transfer within 



the exciting laser pulse (Mairesse et al. , 2010). Another recent implementation 



of this spectroscopy method allowed insights into multi-electron processes in 
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xenon atoms (Shiner et al. 



2011). 



c) Attosecond transient absorption for observing bound- electron wavepackets 
While the photoelectron/-ion and HHG spectroscopy methods discussed above 
crucially require continuum electrons, either directly as observables or as an 
intermediate state to act as a probe interfering with the bound states to be 
analyzed, a third spectroscopy approach — transient-absorption spectroscopy 
(TAS) — is not dependent on free electrons and particularly applicable for the 
spectroscopy of bound-bound transitions. Instead of measuring electrons, it 
measures the spectrally resolved absorption (e.g. absorption lines but also 
non-resonant absorption) of an attosecond pulse transmitted through a coher- 
ently excited sample (e.g. by an intense laser pulse), at different controlled 
time delays. TAS thus does not rely on intermediate or final continuum elec- 
tronic states for dynamical probing, often causing problems in interpretation 
as these states are susceptible to both intense laser fields (typically part of 
the pump or probe process) as well as the long-range Coulomb field of the 
ions that are produced. Continuum electron motion typically does not sepa- 
rate into a Coulomb-only and laser-only perturbation and thus limits direct 
access to bound-bound dynamics and matrix elements without heavily re- 
lying on theory and the validity of suitable models. Photoelectron detection 
techniques also suffer from electron backgrounds produced by the strong laser 
fields alone or by secondary electrons, which do not carry information about 
the dynamics probed by the combined nonlinear action of both the temporally 
separated strong laser field and the weak attosecond pulse. TAS is already 
an established technique in the femtosecond domain, but was only recently 



transferred to gas-phase applications with attosecond HHG light sources (Loh 



et al. , 2007). Using this method, a long-standing scientific question on the 



coherence of different spin-orbit hole states populated in strong-field ioniza- 



tion of atoms (Rottke et al. 1996) could recently be answered ( Goulielmakis 



et al. 2010). In this study, a gaseous but dense (80 mbar) Kr atom sample 



was strong-field ionized by a an intense (3xl0 14 W/cm 2 ) sub-4-fs optical laser 
pulse, promoting the atoms to the 4p7y 2 and 4pg i ionic states, with one elec- 
tron missing from either one of the two spin-orbit split states of the 4p shell. 
A time-delay controlled attosecond pulse exhibiting photon energies around 
80 eV was then sent through the Kr sample to probe these states by transi- 
tions of electrons out of the two spin-orbit-split 3d shell states into the two 
hole states of the 4p shell, resulting in three transitions which are energeti- 
cally located near 80 eV. Characteristic spectral absorption lines were thus 
observed in dependence of the attosecond pulse time delay. Two out of the 
three transitions end up in the same final state, making it possible to probe the 
coherence of the strong-field populated 4p-hole states. A dynamical oscillation 
of absorption strength was observed with a period of 6.3±0.1 fs, correspond- 
ing to the inverse of the energy splitting of these two states, proving their 
coherent excitation during the strong-field pump process. Furthermore, it was 
possible to quantitatively extract the degree of coherence of the excitation, 
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resulting in a value of 0.63 ± 0.17, relatively close to its maximum of unity 
as a result of the short ionization pulse that was used. Recent implementa- 
tion of light-field synthesis methods, by which the ionization pulse could be 



further compressed (Wirth et al. 2011), allowed the temporal reconstruction 



of the spin-orbit wavepacket in real time relative to the inducing laser field. 
Other successful implementations of transient-absorption spectroscopy for the 
measurement of electron dynamics include the laser-induced soft-x-ray trans- 



parency by coupling doubly-excited electronic states (Loh et al. 



2008) and 



measurement of electron wavepacket interference (Holler et al. , 2011[ ) in singly- 
excited He atoms, and the time-domain observation of autoionizing states after 



inner-shell excitation in Ar atoms (Wang et al. 2010). 



6 Outlook 



After more than two decades of HHG, a profound understanding of the under- 
lying mechanism has been reached and numerous applications have been put 
forward. With respect to the latter, attosecond science has been dominating 
and the notion of visualizing atomic structures in motion has been material- 
ized in most beautiful manner. Still numerous challenges remain: The transfer 
of HHG to the multi-keV and MeV range and consequently the generation of 
zeptosecond pulses represents a challenge. In spite of significant progress, laser- 
induced relativistic atomic dynamics still needs to be optimized with regard 
to larger recollision and recombination yields and consequently stronger HHG 
signals. In addition to visualizing nuclear and high-energy processes another 
enormous challenge is the time-resolved spectroscopy of more complex sys- 
tems like e.g. biomolecules. With alternative facilities like free-electron lasers 
and plasma-based schemes becoming increasingly relevant, atomic HHG light 
sources are likely to remain at least competitive due to their unmatched de- 
gree of spatial and temporal coherence, including the synchronization of the 
generated attosecond pulses to optical fields, and due to to their availability 
in numerous laboratories worldwide at reasonable expense. 
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